Calculations of the spectra of superheavy elements E119 and E120 
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High-precision calculations of the energy levels of the superheavy elements E119 and E120 + 
are presented. Dominating correlation corrections beyond relativistic Hartree-Fock are included to 
all orders in the Coulomb interaction using the Feynman diagram technique and the correlation 
potential method. The Breit interaction and quantum electrodynamics radiative corrections are 
considered. Also, the volume isotope shift is determined. A similar treatment for Cs, Fr, Ba + and 
Ra + is used to gauge the accuracy of the calculations and to refine the ab initio results. 

PACS numbers: 32.10.Hq,31.15.am,31.30.Gs,31.30.jf 



I. INTRODUCTION 

There has been great progress in recent years in the 
synthesis of superheavy elements (nuclear charge Z > 
104). Elements up to Z = 118, excluding Z = 117, have 
been produced (see, e.g., Refs. [l|,|2j]), and very recently 
evidence for naturally-occurring E122 was reported 0]. 

Studies of superheavy elements are largely motivated 
by the predicted "island of stability" which occurs due 
to the stabilizing nuclear shell effects. Different nuclear 
models vary in their predictions of the superheavy shell 
structure (see, e.g., [4|). Experimental investigation of 
superheavy elements enables one to distinguish between 
different models. 

Experimental efforts are underway to measure the 
spectra and chemical properties of superheavy elements 
[5|. A number of theoretical works, from the quantum 
chemistry and atomic physics communities, have been 
devoted to these studies (see references in Q and 0). 

Leading relativistic effects grow as (Za) 2 , where a = 
e /he is the fine structure constant, and they become 
very large in superheavy elements. It has been shown 
that these effects lead to a number of interesting features, 
such as level inversion in the spectra of some elements (s- 
levels are pulled in and screen the Coulomb potential seen 
by higher-orbital waves such as <i-waves, thereby pushing 
them out) @. 

In the present work we perform relativistic calculations 
to determine the spectra of superheavy elements El 19 
and E120+. The isotope 292 E120 is predicted to be dou- 
bly magic in relativistic mean- field nuclear calculations. 



II. METHOD OF CALCULATION 

We perform calculations for Cs, Fr, Ba + , and Ra + to 
help gauge the accuracy of the calculations for El 19 and 
E120 + and as a means to reduce the ab initio errors for 
the spectra of these elements through extrapolation. 

At the first stage of the calculations we use the rel- 
ativistic Hartree-Fock (RHF) method. Calculations are 
performed in the self-consistent potential formed by the 
N— 1 electrons in the core {V N_1 potential). A complete 



set of single-electron orbitals is obtained in this way. The 
orbitals satisfy the equation 

hoipa = eo-00 , (1) 
where ho is the relativistic Hartree-Fock Hamlitonian 



ho = ca ■ p + (/3 — l)mc 2 



Ze 2 



V 



N-l 



(2) 



Here = Vdir + K, X ch is the sum of the direct and 

exchange Hartree-Fock potentials, N is the number of 
electrons, N — 1 is the number of electrons in the closed 
core, and Z is the nuclear charge. 



A. Correlations 

The main challenge in calculations of the spectra of su- 
perheavy elements is accurate treatment of correlations. 
We take into account correlations using the correlation 
potential method @. Here, a correlation potential oper- 
ator E is constructed such that its average value for the 
valence electron coincides with the correlation correction 
to the energy, Se a — (a|E|a) 

When the single-particle orbitals are found in the 
Hartree-Fock potential, the many-body perturbation the- 
ory expansion for X starts in second order in the Coulomb 
interaction. There are direct and exchange contributions 
to the correlation potential. Second-order E is calcu- 
lated via direct summation over a discrete set of single- 
particle orbitals. Rather than working with finite sums 
and integrals over the real spectrum, we use finite sums 
over a pseudo-spectrum. We introduce a cavity of radius 
r = 40 a.u and 40 B-splines are used as a basis for the 
functions. 

The ab initio calculations may be improved by in- 
cluding three dominating higher-order diagrams into 
the second-order correlation potential |9j. These are 
(i) screening of the Coulomb interaction, (ii) the hole- 
particle interaction in the polarization operator, and (iii) 
chaining of the correlation potential E. 

In particular, (i) and (ii) are included into the direct 
diagrams of E using the Feynman diagram technique. 
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For the exchange diagrams we use factors in the second- 
order S to imitate the effects of screening. These fac- 
tors are /„ = 0.72, h = 0.62, / 2 = 0.83, / 3 = 0.89, 
fi = 0.94, = 1; the subscript denotes the multipolar- 
ity of the Coulomb interaction. These factors have been 
estimated from accurate calculations of the higher-order 
corrections. The chaining of the correlation potential (iii) 
is included trivially by adding S to the Hartree-Fock po- 
tential. The energies, with correlations included, are so- 
lutions of the equations for the valence electrons, 



(ho + S)V»o = talpa 



(3) 



Further improvements to the wave functions and en- 
ergies may be made semi-empirically through the use of 
fitting factors / (not to be confused with the Coulomb 
screening factors above) placed before the correlation po- 
tential, i.e., 



(ho + fT,)ijj' a = 



(4) 



Factors used for E119 (E120+) are found by fitting to the 
experimental energies of the lighter electronic analogs Cs 
and Fr (Ba + and Ra + ). The use of fitting factors is con- 
sidered a means of including effects, such as higher-order 
correlations, beyond what is included in the ab initio ap- 
proach. 



B. Breit interaction 



where Vjj is the Uehling potential and V g is the potential 
arising from the magnetic formfactor. The potential cor- 
responding to the electric formfactor is divided into low- 
and high-frequency parts, respectively: 



Vi(r) 



B(Z) 



Z i o 5 mc 2 e- Zr/aB 



(8) 



and 



V f (r) = -A(Z,r)-V(r) 

7T 



In (i 2 - 1) + 41n(— 
Za 



dt 



1 



0.5) 



3 1 
~~2 + ¥ 



-Itrra 



(9) 



V(r) is the nuclear potential, the coefficient A(Z, r) = 
(1.071-1.976x 2 -2.128x 3 +0.169a; 4 )mr/(mr+0.07Z 2 a 2 ), 
where x = (Z — 80)a, and as is the Bohr radius. Eqs. 
(|8l9p were determined semi-empirically by fitting to the 
Lamb shifts of high states of hydrogen- like ions for Z=10- 
110. 

This potential is added to the Hartree-Fock potential, 



V 



N-l 



V 



N-l 



rad 



(10) 



It is included in the self-consistent solution of the core 
Hartree-Fock states. Core relaxation, demonstrated to be 
important for the energies of valence p-states, is therefore 
taken into account. 



We go beyond the treatment of the electron-electron 
interaction in the Coulomb approximation, taking into 
account magnetic and retardation effects through inclu- 
sion of the Breit interaction. We use the following form 
for the Breit operator 

B _ 0.1-0.2 + (oi ■ n)(a 2 ■ n) 

where r = nr, r is the distance between electrons, and a 
is the Dirac matrix. 

In a similar way to the Coulomb interaction, we deter- 
mine the self-consistent Hartree-Fock contribution aris- 
ing from Breit. This is found by solving Eq. ([2]) in the 
potential 

V N - 1 = V c + V B , (6) 

where V c is the Coulomb potential, V B is the Breit po- 
tential. 



C. Lamb shift 

Quantum electrodynamics radiative corrections to the 
energies (Lamb shifts) are accounted for by use of the 
radiative potential introduced in Ref . [To| . This potential 
has the form 

14ad(r) - V v (r) + V g (r) + V f (r) + V t (r) , (7) 



III. RESULTS AND DISCUSSION 

We have calculated removal energies for the low-lying 
states s, Pi/2, and Results for Cs, Fr, and E119 are 
presented in Tableland those for the ions Ba + , Ra + , and 
E120 + are presented in Table HT1 We list results in the 
RHF approximation and those with correlations included 
(with dominant diagrams summed to all orders). The 
ab initio results are listed under the column "£" . In 
the column to the right, the percentage deviation from 
experiment is given in brackets. It is seen for Cs, Fr, 
Ba + , and Ra + that there is excellent agreement with 
experiment, with disagreement on the order of 0.1%. The 
largest disagreements are for 7pi/2 for both Fr (0.5%) and 
Ra+ (0.4%). 

In the column "/c s S" in Table U we list the results for 
calculations for Fr and El 19 with the factor /c B found 
by fitting to the measured energies for Cs. It is clear by 
looking at the results for Fr that in all cases the results are 
significantly improved. The deviations from experiment 
are 0.1% or better. It is the same situation for the results 
for Ra + , as can be seen from Table [TT1 

From the trend in the corrections from Cs to Fr, we ex- 
pect that using fitting factors significantly improves the 
accuracy of calculations for El 19. Because use of the 
fitting factors /c s for Fr calculations leads to such good 
agreement with experiment, the fitting factors /ft differ 
only slightly from /c s - This means that extrapolation 
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TABLE I: Removal energies for states of Cs, Fr, and E119 
(cm -1 ), fas, /ft are factors placed before E, found by fitting 
to the energies of Cs and Fr, respectively. Numbers in brack- 
ets are percentage deviations of the preceding result compared 
to experiment. In the last column for E119 we present results 
of another high-precision calculation, Ref. [l|| . Final values 
in the present work for E119 are in column "/fyE". 



State 


RHF 


E 




/csE 


/ Fr E Exp.7 
Other 6 


6s 


27954 


31467 


(0.2) 






31407 


7s 


12112 


12873 


(0.0) 






12872 


8s 


6793 


7090 


(0.0) 






7090 


6^1/2 


18791 


20295 


(0.3) 






20228 


7f>l/2 


9223 


9662 


(0.2) 






9641 


8P1/2 


5513 


5707 


(0.2) 






5698 


6?>3/2 


18389 


19727 


(0.3) 






19674 


7f>3/2 


9079 


9478 


(0.2) 






9460 


8P3/2 


5446 


5623 


(0.1) 






5615 


7s 


28768 


32931 


(0.2) 


32860 (0.0) 




32849 


8s 


12282 


13116 


(0.1) 


13115 (0.0) 




13109 


9s 


6858 


7177 


(0.0) 


7177 (0.0) 




7178 


7P1/2 


18855 


20708 


(0.5) 


20625 (0.1) 




20612 


Spi/2 


9240 


9762 


(0.3) 


9737 (0.0) 




9736 


9P1/2 


5521 


5747 


(0.3) 


5738 (0.1) 




5731 


7P3/2 


17655 


18970 


(0.2) 


18919 (0.0) 




18925 


8P3/2 


8811 


9206 


(0.2) 


9189 (0.0) 




9191 


9f>3/2 


5319 


5494 


(0.2) 


5487 (0.1) 




5483 


8s 


33554 


38954 




38866 


38852 


38577 


9s 


13194 


14087 




14086 


14079 


14050 


10s 


7208 


7534 




7535 


7536 


7519 


8P1/2 


20126 


23445 




23294 


23272 


22979 


9^1/2 


9654 


10453 




10416 


10415 


10365 


10j3l/2 


5709 


6040 




6027 


6018 


5997 


8?>3/2 


16674 


18102 




18046 


18053 


18007 


9P3/2 


8449 


8883 




8863 


8866 


8855 


10j5 3 /2 


5145 


5340 




5332 


5328 


5320 



a Cs data from Ref. [TJ and Fr data from Ref. [TIL 
'Values for E119 are results of calculations, Ref. jl8j |. 



TABLE II: Removal energies for levels of Ba + , Ra + , and 
E120 + (cm -1 ), /ea, /Ha are factors placed before E, found by 
fitting to the energies of Ba + and Ra + , respectively. Num- 
bers in brackets are percentage deviations of the preceding 
result compared to experiment. Final values for E120 + are in 
column "/ Ra E". 



Atom 


State 


RHF 


E 




/BaE 


/ Ra E Exp." 


Ba+ 


6s 


75340 


80834 


(0.2) 






80686 




7s 


36852 


38344 


(0.0) 






38331 




8s 


22023 


22662 


(0.0) 






22661 




6P1/2 


57266 


60603 


(0.3) 






60425 




7P1/2 


30240 


31346 


(0.2) 






31296 




8P1/2 


18848 


19365 


(0.1) 






19350 




6P3/2 


55873 


58876 


(0.2) 






58734 




7f>3/2 


29699 


30718 


(0.1) 






30675 




8f>3/2 


18580 


19060 


(0.1) 






19050 


Ra+ 


7s 


75899 


82034 


(0.2) 


81871 (0.0) 




81842 




8s 


36861 


38454 


(0.0) 


38440 (0.0) 




38437 




9s 


22005 


22677 


(0.0) 


22675 (0.0) 




22677 




7P1/2 


56878 


60743 


(0.4) 


60535 (0.1) 




60491 




8P1/2 


30053 


31297 


(0.2) 


31241 (0.0) 




31236 




9pi/ 2 


18748 


19322 




19306 








7P3/2 


52906 


55771 


(0.2) 


55634 (0.0) 




55634 




8f>3/2 


28502 


29493 


(0.1) 


29451 (0.0) 




29450 




9f>3/2 


17975 


18445 


(0.1) 


1843d (0.0) 




18432 


E120+ 


8s 


83168 


90145 




89964 


89931 






9s 


38468 


40128 




40113 


40110 






10s 


22673 


23357 




23355 


23357 






8P1/2 


60027 


65430 




65141 


65080 






9P1/2 


31121 


32678 




32609 


32604 






10pi/2 


19253 


19945 




19926 


19926 6 






8P3/2 


49295 


52003 




51873 


51874 






9P3/2 


27028 


27993 




27952 


27951 






10p>3/2 


17223 


17691 




17681 


17678 





"Ref. [HI. 

'Final result corresponds to fitting from Ba + , since there is no 
experimental data for Ra + . 



A. Breit and radiative corrections 



of the spectra for El 19 from Fr gives energies that are 
only slightly different from those found from extrapola- 
tion from Cs. Our final results for E119 are found using 
/Fr, presented in Tabled] In the final column of Table U 
we list for E119 results of another high-precision calcula- 
tion [HI and postpone discussion of this work till Section 

mm 

We see the same pattern for the ions, and our final 
results for E120 + are listed under the column "/R a S" . 



Breit corrections were calculated in the self-consistent 
Breit-Hartree-Fock potential and the results are pre- 
sented in Table IIII1 These numbers should be consid- 
ered only as an indication of the order of magnitude of 
the corrections since the correlated Breit corrections may 
be large. For example, in Ref. [IH it was found for 
Cs that account of correlations changes the sign for 6s 
(from 3.2 cm -1 to —2.6 cm -1 ). In that work it was found 
that there is a small suppression due to correlations for 
6pi/2 (7.5 cm -1 to 7.1 cm -1 ) and for 6733/2 it is significant 
(2.9cm -1 to 0.84cm -1 ). 

Our results for quantum electrodynamics (QED) radia- 
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TABLE III: Corrections to removal energies from account of 
the Breit interaction, n is the principal quantum number of 
the ground state. Units are —cm -1 . 



State 


Cs 


Fr 


E119 


Ba+ 


Ra+ 


E120+ 


ns 


3 


6 


35 


14 


26 


82 


(n + l)s 


1 


2 


8 


5 


10 


24 


(n + 2)s 





1 


3 


2 


5 


11 


npi/2 


7 


14 


34 


29 


52 


112 


(n + l)pi /2 


3 


5 


11 


11 


20 


41 


(n + 2)px/2 


1 


2 


5 


6 


10 


20 


np 3 /2 


3 


4 


5 


12 


18 


19 


(n + l)p 3 /2 


1 


2 


2 


5 


7 


8 


(n + 2)p 3/2 





1 


1 


3 


4 


4 



tive corrections are listed in Table HVl alongside results of 
other calculations. We're not aware of other data for the 
ions. As with the Breit corrections, our results should 
only be considered estimates, to give an idea of the size 
of these corrections. They are calculated at the Hartree- 
Fock level, with correlation corrections excluded. (The 
effect of correlations would be to increase the density of 
the valence electrons at the nucleus, thereby leading to 
larger radiative corrections.) 

Moreover, the radiative potential itself was found by 
fitting to states of hydrogen-like atoms for 10 < Z < 110. 
Due to a lack of data, direct fitting for Z = 119, 120 was 
not possible, and it is not clear how well our radiative 
potential would work in this region. 

In Refs. 

0, 

[H|, ratio methods were used to evaluate 
the self-energies. In the former they were found from the 
ratio Ese (Vvp) df/ Eyp, where (Vvp)df is the Uchling 
potential averaged over Dirac-Fock wave functions for the 
neutral system, and Ese and Eyp are self-energy and 
vacuum polarization (Uehling) corrections to the ener- 
gies in hydrogen-like systems. In the latter, the ratio is 
EsE(VU nuc (r))D F /(VU nuc (r))-H, where U nuc is the nu- 
clear potential and ()h denotes averaging over H-like 
states. In Refs. 16] the Lamb shifts are found em- 
ploying rigorous QED in the field of several different ef- 
fective atomic potentials. The numbers in the tables give 
the ranges in the values for the potentials considered. 
Similarly to the current work, in Ref. [17] an effective lo- 
cal potential, mimicking self-energy QED effects, is added 
to the Dirac-Fock potential. 

We see good agreement for the "lighter" atoms, though 
some disagreement for E119. We already mentioned why 
our results should be considered as order of magnitude 
estimates only for the superheavy elements. 

We note that in our work, unlike in all other works 
mentioned, core relaxation is taken into account. This 
is accomplished by including the radiative potential into 
the self-consistent procedure for the core. While this ef- 
fect is relatively small for Cs s levels, it is significant 
for E119. For E119 8s, the Lamb shift changes from 
85cm -1 to 67cm -1 without and with core relaxation, re- 



TABLE IV: Radiative corrections to removal energies. Units 
are -cm" . 



Atom 


State 


This 


Ref. 


Ref. 


Ref. 


Ref. 


Ref. 








[141 


[151 


[161 

1 1U| 


[171° 


[181 

1 10| 


Cs 


6s 


16 


15.5 


14.9^26.6 


12.7-^23.1 


14.1 


18.0 




7s 


4 










4.2 


Fr 




3fi 
ou 


oo.o 


O I . 1 'U 1 . 1 




40 fi 


^o.o 




So 

Go 


q 










9 Q 


ry 


So 


fi7 


1 11 b 

14:1 


11U 




1 


8^ 9 






xo 










99 fi 




8pi/2 


I 










18.2 




8P3/2 


2 










3.7 


Ba+ 


6s 


37 














7s 


12 












Ra+ 


7s 


77 














8s 


24 












E120+ 


8s 


120 














9s 


32 














8P1/2 


5 














8P3/2 


7 













"Self-energies are given in Ref. 11711 : we have added vacuum po- 
larization contributions from Ref. 15] calculated at the Dirac-Fock 
level. 

6 This number is quoted in their later work Ref. [l5f ] without 
explanation; in the original work Ref. 11411 the value is 211 in the 
same units. 



spectively. For p levels the correction is more dramatic, 
although the size of the effect itself is much smaller. For 
El 19 8pi/2, we find the radiative correction to the bind- 
ing energy without and with core relaxation to be 7cm -1 
and 1cm -1 . The effect of the core relaxation is to repel 
the inner electrons (the Lamb shift decreases the bind- 
ing energy), leading to reduced shielding of the nuclear 
Coulomb field at small distances where the radiative cor- 
rections arc determined. 



B. Comparison with other calculations 

We know of only one other work where high-precision 
calculations have been performed for spectra of the super- 
heavy elements studied in this work. Eliav et al. 18] have 
performed coupled cluster calculations for E119 spectra, 
including both Breit and radiative corrections. The re- 
sults of their ab initio calculations are tabulated along- 
side our final (semi-empirical) values in Table HI We see 
that generally there is agreement on the level ~ 0.1%, 
with larger deviations for 8s (0.7%) and 8pi/ 2 (1.3%). 

We investigated the large deviations for levels 8s and 
8pi/2 by calculating the spectra of Cs, Fr, and E119 with 
correlations calculated in the second order of perturba- 
tion theory (£( 2 \ with no higher-order screening or hole- 
particle interactions taken into account). We used fitting 
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factors to mimic higher-order effects, as was done with 
the full correlation potential E, and compared the results 
to those in Table [U The result for 8s obtained in second 

(2) 

order from the fit to Cs, /c 8 ^Eii9> differs from that in 
all orders, /csEeii9, by 0.3% (more bound) and for &P\/2 
the difference is -0.3%. With fitting to Fr spectra (using 
/ Fr and /Fr, respectively, before E^ 2 ) and E) the differ- 
ence is 0.0% for 8s and -0.3% for 8pi/%. For other states, 
the agreement is -0.1% or better. 

For E120 + , very good agreement for s-levels was ob- 
tained using the two approaches with fitting to Ba + and 
Ra + (0.1% or better). For the p-levels, there are larger 
deviations, the largest being 0.5% for 8pi/ 2 - 

Differences in values for the spectra obtained in the two 
approaches gives an indication of the error from missed 
higher-order effects. An estimate of ~ 0.1% error sup- 
ports the detailed consideration below. 



C. Estimate of the accuracy 

Our final results for the superheavy elements do not 
include either Breit or radiative corrections and are listed 
in Tables Q] and [TTJ The reason is that by using factors 
obtained by fitting to measured spectra, it appears that 
some of the Breit and radiative corrections are included. 

Let us consider states of Fr and Ra + with sizeable 
(~ 0.1%) Breit and radiative corrections. For Fr this is 
7s. At the ab initio level, the deviation from experiment 
is 0.25% (column "E"). With the fitting factor / Cs , the 
deviation is reduced to a tiny 0.03%. Moreover, the (esti- 
mated) contribution from Breit and radiative corrections 
(Tables IIIII and IIV[) is much larger than this deviation, 
being 0.13% of the measured energy. For Ra + , looking 
at energies for the states 7s, 8s, 7pi/ 2 , 8pi/2j it is seen 
that we have the same story: the value obtained from 
fitting is everywhere better than the estimated Breit and 
radiative contributions. This strongly supports the argu- 
ment that the use of empirical fitting factors takes into 
account not only the effects of higher-order correlation 
effects, but also the Breit and radiative corrections to 
some extent. 

The question then becomes: can we expect the same 
accuracy for E119 and E120 + as has been demonstrated 
for Fr and Ra + ? Calculations for E119 and E120 + were 
performed in a similar way as for Fr and Ra + and so we 
expect that extrapolation from the lighter to the heav- 
ier systems follows the same pattern we saw from Cs 
and Ba + to Fr and Ra + . However, when we go to the 
heavier systems, there is some difference. For instance, 
the independence of the relativistic, Breit, and radiative 
corrections for the "light" systems is ~ Z 2 , while for the 
superheavy elements this dependence is stronger [14j . It 
means that extrapolation from Fr to E119, for instance, 
is probably not as good as extrapolation from Cs to Fr. 

For Cs, Fr, Ba + , and Ra + it is seen that the largest 
uncertainty in the ab initio calculations comes from the 
unaccounted correlation corrections, these being larger 



than the estimated Breit and radiative corrections. This 
is the case for the higher states for E119 and E120 + , how- 
ever for the ground state these corrections are about the 
same (0.3% and 0.2%, respectively) . (20j We expect that, 
as with Fr and Ra + , use of the empirical fitting factors 
improves the accuracy of the ab initio calculations, and 
accounts somewhat for Breit and radiative corrections. 
We expect our calculations for the superheavy elements 
to be accurate to ~ 0.1%. 



D. Nuclear dependence: volume isotope shift 

For the low s-levels of E119 and E120 + , we have found 
that there is a significant dependence on the root-mean- 
square nuclear charge radius r rms . Our calculations were 
performed using a two-parameter Fermi distribution for 
the nuclear density. The values presented in the pre- 
vious tables were performed with a half-density radius 
c = 8.0 fm and 10-90% width t = 2.0 fm corresponding 
to a rms charge radius r rms « 6.42 fm. Defining the 
volume isotope shift in terms of r rms , 



SE 
~~E 



Sr r 



(11) 



we have found the following values for k for states 8s and 
9s for E119 and E120+ at the RHF level: 



E119 8s : 


k = 


-0.0243 


(12) 


E119 9s : 


k = 


-0.0115 


(13) 


E120+ 8s : 


k = 


-0.0180 


(14) 


E120+ 9s : 


k = 


-0.00936. 


(15) 



A table of values for 



for nuclei Z = 119 and 



Z = 120 calculated in the nuclear Hartree-Fock-BCS ap- 
proximation can be found in Ref. |19j. The values range 
from around r rms — 6.45 fm to r rms — 6.95 fm for the 
very heavy isotopes. For r rms = 6.90 fm, we obtain at 
the RHF level the value 33495 cm -1 for the removal en- 
ergy for E119 8s. The difference between this value and 



that obtained with 



= 6.42 fm, A = -59 cm- 



is 



comparable to the size of Breit and radiative corrections. 

In principle, measurements of the spectra of different 
isotopes of superheavy elements may be used to get in- 
formation about nuclear structure. 



IV. CONCLUSION 

We have performed ab initio calculations of removal 
energies for the the low-lying s and p levels of the super- 
heavy elements El 19 and E120 + . Semi-empirical fitting 
was used to improve the accuracy of the calculations, ac- 
counting for neglected higher-order correlations as well 
as Breit and radiative corrections. The volume isotope 
shift was studied. The accuracy of our calculations is 
estimated to be on the order of 0.1%. 
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